Problems which are well-posed in a generalized sense with applications to the Einstein equations

In the harmonic description of general relativity, the principle part of Einstein equations reduces to a constrained system of 10 curved space wave equations for the components of the space-time metric. We use the pseudo-differential theory of systems which are well-posed in the generalized sense to establish the well-posedness of constraint preserving boundary conditions for this system when treated in second order differential form. The boundary conditions are of a generalized Sommerfeld type that is benevolent for numerical calculation.Comment: Final version to appear in Classical and Qunatum Gravit

Constraint-preserving Sommerfeld conditions for the harmonic Einstein equations

The principle part of Einstein equations in the harmonic gauge consists of a constrained system of 10 curved space wave equations for the components of the space-time metric. A new formulation of constraint-preserving boundary conditions of the Sommerfeld type for such systems has recently been proposed. We implement these boundary conditions in a nonlinear 3D evolution code and test their accuracy.Comment: 16 pages, 17 figures, submitted to Phys. Rev.

Modeling the Black Hole Excision Problem

We analyze the excision strategy for simulating black holes. The problem is modeled by the propagation of quasi-linear waves in a 1-dimensional spatial region with timelike outer boundary, spacelike inner boundary and a horizon in between. Proofs of well-posed evolution and boundary algorithms for a second differential order treatment of the system are given for the separate pieces underlying the finite difference problem. These are implemented in a numerical code which gives accurate long term simulations of the quasi-linear excision problem. Excitation of long wavelength exponential modes, which are latent in the problem, are suppressed using conservation laws for the discretized system. The techniques are designed to apply directly to recent codes for the Einstein equations based upon the harmonic formulation.Comment: 21 pages, 14 postscript figures, minor contents updat

Testing the well-posedness of characteristic evolution of scalar waves

Recent results have revealed a critical way in which lower order terms affect the well-posedness of the characteristic initial value problem for the scalar wave equation. The proper choice of such terms can make the Cauchy problem for scalar waves well posed even on a background spacetime with closed lightlike curves. These results provide new guidance for developing stable characteristic evolution algorithms. In this regard, we present here the finite difference version of these recent results and implement them in a stable evolution code. We describe test results which validate the code and exhibit some of the interesting features due to the lower order terms.Comment: 22 pages, 15 figures Submitted to CQ

Well-posed initial-boundary value problem for the harmonic Einstein equations using energy estimates

In recent work, we used pseudo-differential theory to establish conditions that the initial-boundary value problem for second order systems of wave equations be strongly well-posed in a generalized sense. The applications included the harmonic version of the Einstein equations. Here we show that these results can also be obtained via standard energy estimates, thus establishing strong well-posedness of the harmonic Einstein problem in the classical sense.Comment: More explanatory material and title, as will appear in the published article in Classical and Quantum Gravit

Boundary conditions for coupled quasilinear wave equations with application to isolated systems

We consider the initial-boundary value problem for systems of quasilinear wave equations on domains of the form $[0,T] \times \Sigma$, where $\Sigma$ is a compact manifold with smooth boundaries $\partial\Sigma$. By using an appropriate reduction to a first order symmetric hyperbolic system with maximal dissipative boundary conditions, well posedness of such problems is established for a large class of boundary conditions on $\partial\Sigma$. We show that our class of boundary conditions is sufficiently general to allow for a well posed formulation for different wave problems in the presence of constraints and artificial, nonreflecting boundaries, including Maxwell's equations in the Lorentz gauge and Einstein's gravitational equations in harmonic coordinates. Our results should also be useful for obtaining stable finite-difference discretizations for such problems.Comment: 22 pages, no figure

Finite difference schemes for second order systems describing black holes

In the harmonic description of general relativity, the principle part of Einstein's equations reduces to 10 curved space wave equations for the componenets of the space-time metric. We present theorems regarding the stability of several evolution-boundary algorithms for such equations when treated in second order differential form. The theorems apply to a model black hole space-time consisting of a spacelike inner boundary excising the singularity, a timelike outer boundary and a horizon in between. These algorithms are implemented as stable, convergent numerical codes and their performance is compared in a 2-dimensional excision problem.Comment: 19 pages, 9 figure

Hyperbolicity and Constrained Evolution in Linearized Gravity

Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint equations remain solved under the action of the evolution, and one approach is to simply monitor them ({\it unconstrained} evolution). Since computational solution of differential equations introduces almost inevitable errors, it is clearly "more correct" to introduce a scheme which actively maintains the constraints by solution ({\it constrained} evolution). This has shown promise in computational settings, but the analysis of the resulting mixed elliptic hyperbolic method has not been completely carried out. We present such an analysis for one method of constrained evolution, applied to a simple vacuum system, linearized gravitational waves. We begin with a study of the hyperbolicity of the unconstrained Einstein equations. (Because the study of hyperbolicity deals only with the highest derivative order in the equations, linearization loses no essential details.) We then give explicit analytical construction of the effect of initial data setting and constrained evolution for linearized gravitational waves. While this is clearly a toy model with regard to constrained evolution, certain interesting features are found which have relevance to the full nonlinear Einstein equations.Comment: 18 page

The Well-posedness of the Null-Timelike Boundary Problem for Quasilinear Waves

The null-timelike initial-boundary value problem for a hyperbolic system of equations consists of the evolution of data given on an initial characteristic surface and on a timelike worldtube to produce a solution in the exterior of the worldtube. We establish the well-posedness of this problem for the evolution of a quasilinear scalar wave by means of energy estimates. The treatment is given in characteristic coordinates and thus provides a guide for developing stable finite difference algorithms. A new technique underlying the approach has potential application to other characteristic initial-boundary value problems.Comment: Version to appear in Class. Quantum Gra